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The  neoclassical  transport  coefficients  in  the  low  collision  frequency  regime  have  been 
investigated  by  many  authors  (HINTON  and  HAZELTINE,1976,  HIRSHMAN  and  SIG- 
MAR.1981).  Using  the  large  aspect  ratio  expansion,  ROSENBLUTH  et  al.(1972)  obtained 
the  transport  coefficients  from  a  variational  principle.  The  aspect  ratios  for  the  present  and 
the  planned  future  thermonuclear  reactors  are  not  so  large  that  the  finite  aspect  ratio  correc¬ 
tion  is  important.  Recently,  BOLTON  and  WARE  (1983)  calculated  the  ion  thermal  conduc¬ 
tivity  without  assuming  the  large  aspect  ratio.  They  solved  the  drift  kinetic  equation  numeri¬ 
cally  by  using  the  model  collision  operator  derived  by  HIRSHMAN  and  SIGMAR  (1976). 
The  finite  aspect  ratio  correction  is  shown  to  increase  the  ion  thermal  conductivity  by  a  factor 
of  two  in  the  intermidiate  aspect  ratio  case.  In  this  paper,  wc  introduce  a  new  approximate 
collision  operator.  Using  this  collision  operator,  we  calculate  the  ion  thermal  conductivity 
and  the  ion  distribution  function  in  the  banana  regime  for  an  axisymmetric  toroidal  plasma  of 
arbitrary  aspect  ratio. 

Let  us  introduce  the  flux  coordinates  (  tjf.'f.x)  where  ij»  is  the  poloidal  flux,  ^  the  toroidal 
angle, and  \  the  poloidal  angle-like  variable.  We  consider  a  plasma  composed  of  electrons 
and  a  single  kind  of  ions  in  an  axisymmetric  magnetic  field 

B  =  /V<j+  V^xVij/,  (1) 

where  I  is  a  function  of  i|>  only.  The  mass  and  the  charge  of  an  electron  are  denoted  by  mr 

and  ~  e,  and  those  of  an  ion  by  m,  and  et. 

The  basic  equation  to  describe  the  neoclassical  transport  is  a  drift  kinetic  equation  for  a 
guiding-center  distribution  function.  In  many  situation  of  practical  interest,  the  guiding- 
center  distribution  function  fi  for  ions  is  slightly  perturbed  from  the  local  Maxwellian  distri¬ 
bution  function  /i0.  We  assume  that  the  perturbation  fn  due  to  the  ion-temperature  gradient 
frora/l0  is  the  first  order  quantity  in  a  parameter  b  =  pip/L,  where  p,p  is  the  poloidal  Larmor 
radius  of  an  ion  and  I.  is  the  plasma  scale  length.  Then,  the  distribution  function/),  satisfies 
the  linearized  drift  kinetic  equation 

\  1b~^x  "  CM  =  “  5/2)^/.o> 


(2) 


where  we  have  chosen  the  energy,  the  magnetic  moment  p.=  v*!2B  and  ct= vir/|v„|  as  indepen¬ 
dent  variables  instead  of  the  velocity.  Here.v,  =  vB/|B|,jc=  v/v, , vi  =  y/2 T/in t ,T\= d T/d i(/  ,TI 
the  temperature  of  ions,  =  (Im^/e^-V^/B)  the  radial  component  of  drift  velocity 

and  Cuthe  linearized  Fokker-Planck  operator  for  collisions  between  ions.  We  have  neglected 
collisions  with  electrons  in  the  above  equation  (2)  since  the  mass  ratio  m,/ml  is  very  small. 

In  the  banana  regime,  we  expand  the  distribution  function  as  fnss/l0)  +  vuf}l)+  •  •  •  > 
where  vji=4irnie/lnA/m?v<3  is  the  ion-ion  collision  frequency.  Then,  the  function  /J0)  and/fl) 
satisfy  the  equations 

v*  "  V«-*K**  “  (3) 

v*  lf"7x~  "  C“^#>)  ■  °-  (4) 

The  solution  to  the  eq.(3)  is  given  by 


yjo,  _  _  f^.(j4  _  5/2)fn  +  G(x,X,o)|  (5) 

where  </-»{, /v,\=(l-<72yfl.  The  function  G  is  obtained  from  the  following  solvability  condi¬ 
tion  for  the  eq.(4): 

S-*irC»(lt(xl-5/2V»  +  G)  -  0.  (6) 

vn  °x  ° 

Here  the  integration  means 

[or  k<k-  (7) 

where  'Kc^VBmtx,  and  Xi(M  and  XjM  are  bounce  points  of  trapped  particle.  As  is  well 
known  (ROSENBLUTH  et  al.,1972),  the  function  G  becomes  identically  zero  in  the  trapped 
particle  region. 

The  heat  flux  across  a  magnetic  surface  is  given  by 


Qt  -  <Jdv|f»,v’  - 
-  -2— L<  fdvf(x>  - 1 ■K.m  >. 


(8) 
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x-4-T  [C'uixtf-SWJ  +  (C*  +  vtt(*))*(*)], 

vuW 

where  8(Xf-X)  is  the  Heaviside  function.  The  function  K(x)  satisfies  the  equation 
vu(x)K(x)  -  /e(C2  +  va(*))jr(*)  -  f'CttxtfS/lVn), 

where 


/  „  l<fl2>fc  XdX - 

4  <Vl  -XB> 

This  equation  is  derived  from  (15)  and  (16).  Using  the  expression  (16),  we  obtain 


The  heat  flux  is  also  written  in  terms  of  K(x)  as 


(16) 

(17) 

(18) 

(19) 


L,  =  -<—><B2>fg,‘dxx3(x2  -  5/2 )C/((x(x2  -  5/2)e',J) 

+  ^^-fJdxx3(S/2  -  x2)Cl)(K(x)e~*2),  (21) 

where  X(x)  =  X(x)//l0,C}l=Cll/vu,  and  T,  =  3V2it/4v(j  . 

We  calculate  the  function  K(x)  by  means  of  the  Chapman-Enskog  method.  Let  us  expand 

ff  _  j 

the  function  K{x)  as  K{x)  =  x  £  (**)/«  .where  is  the  Sonine  polynomial  of  order  3/2. 

j-o 

By  substituting  this  expansion  into  (17),  we  obtain  the  following  equation  for  a,: 

-  V  [(!  -  fc)<JFuW>  -  fc<j&\k>]at  -  -  fc<j\Cu\l>,  j  =  0.1.2...JV-U22) 

»-o  *- 

where  vu{x)ssvu{x)/vu  .  Here 


<;|Cl/’>  =  /0  dxx3S%\ (x2) C (xS$ (xJ)e - *2) ,  (23) 

where  C  =  C,‘  or  v„(x).  The  matrix  elements  <j\C\j‘>  are  calculated  by  using  the  generating 

functions  (KANEKO.1960,  TAGUCHI, 1982,1983). 

In  the  limits  of  unit  aspect  ratio  and  large  aspect  ratio,  the  ion  thermal  conductivity  can  be 
obtained  analytically.  The  thermal  conductivity  for  the  unit  aspect  ratio  (  /,  =  0  )  becomes 
/.,  =  <B>2<l/flJ>/V2  since  K(x)= 0.  This  result  agrees  with  the  one  of  HAZI-LTINF.  et 


al,(1973).  The  function  K(x)  for  the  uniform  magnetic  field  (/c  =  l  )  is  given  by 


K{x)  =  |*(f  -  * 2)  ~ 

By  means  of  the  relation 


<0lv,vWll> 

<o!v„W|o> 


fo'd xe-*2x>(S/2  -  x*)Cl(K)  =  <l\Cj,\l> 

+  (1  -fc  y/cf0  dx e-JV(5/2  -  (25) 

we  can  obtain  the  thermal  conductivity  in  the  large  aspect  ratio 

(<f>  -  7k) +  (1  -  s  (1  -  «■  •  (26) 

“  ~  ~  +  <lp,W|1>  S  °  326,  (27) 
This  agrees  with  the  result  of  ROSENBLUTH  et  al.(1972). 

Next,  we  derive  the  simple  expression  for  the  ion  thermal  conductivity.  The  lowest 
approximation  for  £(x)  is  given  by  N=2.  Then.the  approximate  thermal  conductivity 


becomes 


M„>  _ l _  +  4_<8>L,  ll/5- 

'  V2  B 2  <B2>  v2  <b2>  1  +  p/c 


P  =  "1+  vSTs  L17-  (29) 

As  will  be  shown  later,  this  is  in  very  close  agreement  with  the  more  rigorous  one  obtained 
by  N=50.  Note  that  the  approximate  thermal  conductivity  L\ap)  for  the  limits  of  unit  aspect 
ratio  and  large  aspect  ratio  agrees  with  the  exact  one  L,. 

Finally,  we  practically  calculate  the  thermal  conductivity  L,  and  the  ion  distribution  func¬ 
tion  /fu)  in  a  typical  magnetic  field  of  tokamak  with  circular  cross  section  ,  i.e.  , 


B  =  (B„BvBt)  =  (0  ,ByJh,Bjh),  (30) 

h  -  1  +  ecosfl ,  t  =  r/R,  (31) 

Here  r  and  0  are  the  polar  coordinates  in  a  plane  perpendicular  to  the  magnetic  axis  and  R  is 

the  major  radius.  For  this  magnetic  field,  fr  is  computed  numerically  and  is  shown  in  Fig.l. 

The  dashed  curve  represents  the  analytic  expression  /.  =  1-  1.035V27  for  small  e.  In  Fig. 2, 

we  depict  the  thermal  conductivities  L,  =  VlL/Vt  and  Ftf  =  \^2l as  a  function  of  e. 
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From  this  figure,  the  approximate  thermal  conductivity  L°p  is  shown  to  agree  well  with  the 
more  rigorous  one  Lr  CHANG  and  HINTON  (1982)  derived  the  following  approximate 
expression, 

L{  =  (0.66  +  1.88 Vt  -  \.54t)<B>1<UB1>  (32) 

This  expression  was  obtained  by  using  the  BOLTON  and  WARE’S  result  for  e  =0.3,  and  the 

results  of  ROSENBLUTH  et  al.  for  the  large  aspect  ratio  and  of  HAZELTINE  et  al.  for 

the  unit  aspect  ratio.  For  comparision,  we  also  plot  the  values  of  this  expression  in  Fig.2  by 

dashed  curve.  The  ion  distribution  function  /^0)  is  easily  computed  from  (19).  In  Figs. 3(a)- 

(d),  we  shows  the  contours  of  /<0)//l0  for  0  =0,  and  (a)  e  =0.1,  (b)  €  =0.3,  (c)  e  =0.5.  The 

contour  levels  are  evenly  spaced.  The  dashed  lines  represent  the  boundaries  between  the 

trapped  and  the  untrapped  region.  In  Fig. 4,  the  distribution  functions 

/J0> = /J0)//J0(— mic/v47*/Ve<7*i)" 1  are  plotted  as  a  function  of  v(/v,  for  0  =0,  =0,  and  e 

=0,0.1, 0.3  and  0.5. 

In  conclusion,  we  have  presented  the  new  approximate  collision  operator.  Using  this 
operator,  we  have  formulated  the  method  for  calculating  the  ion  thermal  conductivity  and  the 
ion  distribution  function  in  the  banana  regime.  Our  method  can  be  used  for  the  axisymmetric 
toroidal  plasma  of  arbitrary  aspect  ratio. 
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F1GURE  CAPTIONS 

Fig.l  Values  of  fc  as  a  function  of  t. 

Fig.2  The  ion  thermal  conductivities  L,  and  L‘r  as  a  function  of  €.  The  dashed  curve 
shows  the  expression  (32)  obtained  by  CHANG  and  HINTON. 

Fig.3  The  contours  of  /J°V/j0  for  0  -0  and  the  inverse  aspect  ratios  (a)  6  =0.1,  (b)  e 
=0.3,  (c)  c  =0.5.  The  contour  levels  are  evenly  spaced.  The  dashed  lines  represent 
the  boundaries  between  the  trapped  and  the  untrapped  region. 

Fig.4  The  ion  distribution  functions  fi0)  are  plotted  as  a  function  of  v;,  for  0  =0,  ^  =0, 
and  t  =  0,0. 1,0.3  and  j.5. 


